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A path integration formulation for the finite density and temperature problems is shown to be 
consistent with the thermodynamics using an 8 component "real" representation for the fermion 
fields by applying it to a free fermion system. A relativistic quantum field theory is shown to be 
smoothly approached at zero temperature by a real-time thermal field theory so derived even at a 
finite density. The analysis leads to a new representation for the fermion fields which is shown to 
be inequivalent to the conventional 4 component theory at the quantum level by having a mirror 
universe with observable effects and to be better behaved at short distances. 
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. The relativistic quantum field theory (QFT) at finite density and temperature is currently a viable tool to study 
the real-time processes of a relativistic system under extremely conditions like in heavy ion collisions, in astronomical 
and cosmological processes, etc. A sufficiently general formulation of the relativistic QFT at conditions different from 
zero temperature and density is necessary so that it can handle the cases of non-perturbative spontaneous symmetry 
breaking and particle production and is consistent with thermodynamics in real-time. Such an extension is not trivial. 
One of the non-perturbative treatments of the relativistic QFT is based upon Feynman-Matthews-Salam (FMS) 
\ path integration jjj. The FMS formalism expresses time evolution between initial and final states in terms of an 
^ ■ integration over paths that connect the initial and final states with a weight determined by the action of each particular 
' path. The results of the path integration at a formal symbolic level are not uniquely defined in the Minkowski space- 
, time due to the presence of singularities. Their uniqueness is determined by the fact that the derivation of the FMS 
path integration representation of the evolution operator implies a particular ordering of the intermediate states which 
defines how the singularities of the formal results are going to be handled. Thus the FMS path integration formalism 
\q 1 is defined by not only the formal expressions that contain singularities but also by the "causal structure" following 

| the time ordering of the physical intermediate states. 
""^ i One form of the thermal field theory (TFT) was developed in Ref. in the Euclidean space-time with the time 
variable playing the role of the inverse temperature. Such a formalism, which is based on the periodic (antiperiodic) 
JL , boundary condition for bosons (fermions) fields, provides the correct thermodynamics for a noninteracting system. 
The periodic (antiperiodic) boundary condition for the field operators is explicitly written as 

i>(t - i/3) = T)e-^m, (1) 



where (3 is the inverse temperature, 77 = ±1 for bosons and fermions respectively and /1 is the chemical potential. 
The physical responses of the system to external stimulations that happen in real time can be obtained by a proper 
analytical continuation. Another formulation of the same problem based upon Eq. [I] can be derived by a distortion 
of the Matsubara time contour, which goes straightly from to to a contour that contains the real time axis 
extending from negative infinity to positive infinity and returning contour below the real time axis somewhere between 
and — i(3 that parallels the real time axis (see, e.g., Ref. ||). 

Consistency requires that the results of the TFT go smoothly to that of a relativistic QFT in real time at zero 
temperature (or (3 — > 00). The question is whether or not it actually happens, especially for fermions. To investigate 
such a question, let us study the thermodynamics of a free fermion system at finite density and zero temperature 
using relativistic QFT defined by a FMS path integration formalism. 

The thermodynamics of a system with variable fermion number density is determined by the grand potential f2 
defined in the following 

e -W = K e -0(H o -»N) > (2) 

where N is the fermion number density operator and Hq is the (free) Hamiltonian of the system. The result for f2 is 
well known from elementary statistical mechanics; it is U — TS — fj,N with U the internal energy, T the temperature, 
S the entropy and N the average particle number of the system. 

In the quantum field theoretical investigation of the same system at zero temperature, there is another representation 
for the right hand side (r.h.s.) of Eq. || in terms of Feynman-Matthews-Salam |lj path integration 
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lim < e- f3K >= const x / [Dip][Dip]e- SE , (3) 

where K = Hq — uN, "const" is so chosen that 17 = at zero fermion density. Se is the Euclidean action of the 
system; extending the rules of Ref. Q to finite density cases, it can be obtained from the Minkowski action for the 
system by a continuous change of the metric. The result is 

„5 



d A xip (ip + fij 5 -m)ip, (4) 

where a; M is the (4-dim) Euclidean space-time coordinates, 7^ = i-f 5 and m is the mass of the fermion. ft is obtained by 
identifying f3 with X4 (the Euclidean time) that confines the system in the temporal direction before the thermodynamic 
limit. From Eq. ^, the path integration can be carried out immediately. The result |^| is 

ft' = — ^lim — LnDet7^. (ift + /17 5 — to) — ln(const) 

v^co J 2ir^ 2vr p 2 + to 2 V ; 

with the Euclidean 4 momenta p^ = (pi,p) and p^_ = (p 4 — ifi,p). The result for ft' is finite despite the infinite 
integration in momentum. For example, in case of zero mass PV = —ft' = V/i 4 /I2ir 2 , which is expected from the 
elementary statistical mechanics. 

The Euclidean energy pi can be treated as a complex number so that the pi integration in the complex Minkowski 
energy po plane is represented by contour Ce in Fig. I]. In the study of the real-time processes of the system, Minkowski 
action should be used. The original po integration contour Cfms corresponding to the FMS path integration in the 
Minkowski space-time is also shown in Fig. 0. There are also a pair of contours C+ and C_ shown in Fig. ||, called the 
quasiparticle contour, that contribute to the real-time response of the system in the (3 — > 00 limit. Contours Cfms, 
Ce, C + and C_ in Fig. || belong to the same topological class of contours having the same FMS causal structure. 
Consistency requires the equivalence of the set of contours Cfms, Ce, C+ and C_ for the physical quantities. Eq. |5| 
unfortunately fails to meet this requirement due to the fact that the imaginary part of the logarithmic function falls 
off as O(n/\po\) on the physical po sheet (the imaginary part of the logarithmic function on the edges of its branch 
cuts on the physical sheet is shown in Fig. [I]). This causes the results obtained from doing the po integration on the 
above mentioned set of contours different from each other since the integration on the sections of the large circle of 
contours that connecting them has a non- vanishing value. 

Explicit computation of Eq. ^ on contour C+ where only the imaginary part of the integrand with value Hn con- 
tributes (see Fig. ^), reduces to a form differs quite drastically from the expected grand potential of thermodynamics, 
which is ft = U — fiN. In fact, it diverges. 

The cause of the above mentioned problems is found to be related to the asymmetric nature of the 4 component 
representation of the fermion fields with respect to particles and antiparticles. An 8 component "real" representation 
for the fermion fields can be used ||. For a spin 1/2 particle with one flavor, its 8 component spinor ^ can be written 
as \I/ = (ipi,ip2 ) T with the reality condition given by 

* = 7W* T , (6) 

where superscript "T" represents the transpose. Here M. = 0\ ® C with C the ordinary charge conjugation operator 
and 0\ the first of three Pauli matrices 0^2,3 acting on the upper (ipi) and lower (-02) 4 components of 'J. In the 
Minkowski space-time, \& transforms in the same way as ipi or -02 under the Poincare group, which commutes with Oi 
(i = 1, 2, 3). Covariant operators can thus be constructed in the same way as in the 4 dimensional representation. For 
example, the bilinear terms in \1/ can be written as ^T®0\^, where T is any covariant matrix in the 4 dimensional Dirac 
space, A = {0, 1, 2, 3} and Oq = 1. Only the matrices satisfying the antisymmetric condition (M.T®0\) T — —M.T®0\ 
can be selected since $ is a collection of anticommuting Grassmanian numbers in the path integration language. 
In terms of the invariant action Se (Eq. |]) for a free massive fermion system in the Euclidean space is 

Se = § / d A x^> (ift + uj 5 3 - to) * and Eq. | becomes lim^oo < e~ pK >= const x J[D^]e^ SE where the ^ 
degrees of freedom are not functionally integrated due to the reality condition Eq. ^|. The path integration can again 
be easily evaluated: 

J (2tt) 4 (p 2 + m 2 ) 2 ' 
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with the same order of integration as Eq. ||. Here pt = (P4 + i(-i,p)- Eq. |?] and Eq. || have an identical value on 
the contour Ce- They differ on contours C± because the large po behavior of the imaginary part of the logarithmic 
function in Eq. [7] is of order 0(/j, 2 /|j?o| 2 ) on the physical sheet, which allows the equivalence between the set of 
contours Cfms, Ce, C+ and C_. By following the C+ contour shown in Fig. H it is simple to show that the resulting 
r.h.s. of Eq. ^ is finite and unique, namely, U — fiN. It is the zero temperature grand potential for a free fermion 
system at density n — N/V expected from thermodynamics. 

The correct grand potential for a free system in real-time can be obtained using a 4 component representation for 
the fermion fields in the TFT. The standard method do a Matsubara summation over discrete (imaginary) energies 
implied by Eq. [l| for the 4 component Dirac spinor ip. Such a boundary condition for ip results, however, in different 
analytic structure for the path integration than that for the FMS path integration in the real time. The grand 
potential in the Matsubara formalism is 

n = - hm 2 Z f d3p 
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v^oo /3 J (2tt) 

~ m 2 +p 2 -[i(2n+l)nP^+n} 2 

> In 5 — . (8) 

m 2 +p 2 - [i(2n + l)^- 1 ] 2 

The sum over Matsubara frequencies can be evaluated by using a contour integration. For this case, it is 
2itif3~ 1 J2 n f( z n) — f c dztaah\j3z/2]f(z), where f[z n = i(2n + l)7r/3 _1 ] is the logarithmic function in Eq. ||, Co 
is shown in Fig. || and the proper f(z) that gives the correct thermodynamics is f(z) = ln[m 2 + p 2 — (z + fj,) 2 ] — 
ln[?7i 2 + p 2 — z 2 }. Since f(z) approaches to zero fast enough in the \z\ — ► 00 limit and is analytic in the complex z 
plane excluding the real axis, the integration over z along contour Co is equivalent to the sum of integration over 
contours C+ and the negative of C_, which gives Q = U — TS — fiN = —PV expected from thermodynamics. 

The reason why the FMS path integration evaluation for the r.h.s. of Eq. || using a 4 component fermion fields 
fails to give the correct thermodynamics at zero temperature can be understood by a comparison of Figs. |l| and |^. 
The causal structure of the FMS path integration restricts the po integration for the partition function to be within 
the class of contours shown in Fig. |l|. It lacks the contour pair C+ and C_ implied by the boundary condition Eq. |l| 
shown in Fig. which are necessary components to given the correct thermodynamics. With the 8 component "real" 
spinor 'J to represent the fermion fields, the effects of the C_ are provided by the lower 4 component 1/12 of f . 

We are ready to develop a real-time field theory at finite density and temperature using the 8 component "real" 
representation for the fermion fields. The boundary condition Eq. |l| is written in an equivalent form. For the fermions 

interested in this study, it can be expressed as ^(t) = —e^ 3Kl i'(t)e~^ K . It is equivalent to Eq. |l| due to the fact that 
the conserved N commutes with the total Hamiltonian H, which allows the factorization of the action of f3fiN and 
(3H in the exponential of (3K with the exp(/3/x) factoiQ in Eq. |l| the result of the action of exp(/3^iV) and exp(—f][iN) 
on if) from left and right. To handle non-perturbative symmetry breaking cases, the commutativity of H and N shall 
not be imposed at this level of development but rather at later stages as dynamical constraints. 

The real-time thermal field dynamics in the 8 component "real" representation for fermion fields can be constructed 
along the same line as that of Refs. g), Q, 0- ||. Compared to those developments, the only formal differences of 
the present approach are 1) the Kubo-Martin-Schwinger (KMS) boundary condition for the contour propagator do 
not contain the exp(/3/z) factor in the present approach, 2) the effects of \i is hidden in the energy variable within the 
propagator (namely, the time evolution is generated by K not H) here and 3) the symmetry factor for a Feynman 
diagram in a perturbation expansion is different from the 4 component representation for the fermion fields due to 
Eq. ^, which also allows the Wick contraction between two ^s ||. The matrix form of contour propagator^ for a 
fermion is 

S{p) = M ( 5 Q + 5 °_ \ M, (9) 

where the retarded and advanced propagators S + and S~ are S ± (p) = i/(i> + 7VO3 - m ± i7°e) and M = 
sgn(p )\/l n (Po) — i^y/nipo) with n(e) = l/[exp((3e) + 1] and 02 the second Pauli matrix. The poles of the 
above propagator at zero temperature are located above the contour Cfms if they are negative and below it oth- 
erwise. It can be compared to the conventional real-time approach with the correct causal structure pO| , pl] | at zero 



1 It should be exp(/3/iOa) in the new representation. 



2 The returning time contour parallels the real time axis is chosen to be located half way between and — ip. 
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temperature, in which the poles of the fermion propagator are located above the po integration contour if they are 
on the left of p, or else below it. The real-time partition function obtained here for the system is identical to the 
Matsubara method since it has a better large po behavior at finite density. 

Let us turn to the study of local bilinear operators constructed from two fermion fields of the form O = ip{x)Tip[x) 
in a finite density situation where Y is certain matrix acting on ip. A product of two field operators is in general 
singular and non-unique. A definition of such a potentially divergent product is O = lim^^o tp(x + 5)Tif)(x), where 
6^ is a 4-vector with, e.g., S < 0. The ground state (vacuum) expectation value of O can thus be computed from 
the fermion propagator S by closing the integration contour on real po axis in the lower half of the complex po plane. 
Since there is no poles and cuts off the real po axis except the ones on the imaginary po axis due to the thermal 
factor, which should be excluded (see Fig. ||), it can be deformed to Cr to include the poles and cuts of the retarded 
propagator S + : 

with "tr" denoting the trace over internal indices of the fermion fields. Cr reduces to C+ as (3 — > 0. For a free system, 
there are only poles in S 11 so that Eq. 10 can be expressed as 



k d > = i S ~ ^°)] trRes5+ (^) r > ( n ) 

where the sum is over all poles pi of S + (p) and "Res" denotes the corresponding residue. The difference between the 
8 component theory here and the 4 component one also manifests here. For example, the conserved fermion number 
density n corresponding to current j» = i*7 M 3 * is n = ^ dpp 2 - /^ +) with /^ ±} = \/[ e P( E P±n) + i], it 

is same as the one in elementary statistical mechanics. On the other hand, n computed in the conventional approach 
using such a point split definition of fermion number density is divergent and different from each other for S > and 
5q < 0, which means that it is not even consistent with relativity it starts with for a space like 6^ whose time component 
can has different sign in different frames. It can be made finite and unique only after an arbitrary subtraction. 

The particle content of the 8 component theory can be found by studying pole structure of the fermion propagator 
together with the consideration of the FMS causal condition. The time dependence of the propagator at zero tem- 
perature is S(t,p) = f CpMS ^' e ~ Wot S{po 1 p), which, in case of t > 0, can be evaluated on the contour C + with a 
result 

S(t, p) = 0{E P - M )Ai+e-^-^* + OQi - E p )K 2 -e- l ^- E ^ + A^-^^)* (12) 

and, if t < 0, can be evaluated on the contour C- to obtain 

S(t,p) = K]- e ^ + ^ t + 9<ji - E^Al+e 1 ^-^' + 0{E p - ^A^e^^-^*. (13) 

Here A£± = P r (±j°E p -7 p + m)/2E p , Pi = (1 + 3 )/2 and P 2 = (1 - 3 )/2 are projection operators. 

The FMS causal structure (or boundary condition) in the present theory can be simply putted as: 1) excitations 
with po > that propagate forward in time correspond to particles or antiholes and 2) those with po < that 
propagate backward in time correspond to antiparticles or holes. The quantization of the 8 component fermionic field 
then follows naturally. In case of zero density (// = 0), ^(x) can be written as 



= VfF \u rps e- ip - x b rps + V rps e tp - x dK P s] ■ (14) 



2E„ 

rps ^ 



Here s is the spin index, r = 1,2, U rps and V rps are 8 component spinors that satisfy {j> — m)U rps — 0, O^U rps = 
(3 — 2r)U rps , (j/> + m)V rps — 0, O^V rps = (3 — 2r)V rps and the non vanishing anticommutators between b rps and d rps 
arc |S rps , b\., p , s , I = 2E p 5 pp >5 ss >5rr' , |drps, d\., p , s , j = 2E p S pp iS ss ' <5rr', which lead to a canonical quantization of <!». 
The particle states are obtained from the vacuum state | 0), which is annihilated by the bs and ds above, by the action 

of <f and ^ on it. Albeit proper counting has already been given to the intermediate states during its evolution [||, 
these particle states are not all detectable since the reality condition Eq. |^ has to be imposed on the detector. The 
operator reality condition for the fermion field 'J are given by the mirror reflection operator M, which is defined by 
M^> T M = M<& T with M 2 = I, where <& is obtained from ^> by the following replacement: Mb rps M — d rps . The 



4 



operator reality condition corresponding to Eq. || is then * = M<& T . The vacuum state satisfies P p h ys \ 0) =| 0). So 
the physical states satisfying the operator reality condition are those ones symmetric in r = 1 and r = 2 projected 
out by Pphys = (I + M)/2. In finite density case, the annihilation and creation operators for the holes (in the Fermi 
sea) has to be introduced, which is straightforward. 

Compared to the 4 component theory, which has only one fermionic excitation excitation mode with an effective 
energy e_ = E p — [i {p > 0) for particle, which is the one with r = 1 in Eq. [l^, Eq. [l^ has a new fermionic excitation 
mode with r = 2 and an effective energy e + = E p + [i possessing a common 3 momentum range, spin structure and 
opposite charge against the r = 1 one. This excitation mode is much larger in energy than the r = 1 one (e_) in a 
non-relativistic system like the electron gas, it can be comparable in relativistic systems like the quark gluon plasma 
where the bare mass of the fermion is very small relative to its momentum. In addition, there is another particle mode 
with r = 2 and effective energy fj, — E p . Its effects will be further discussed in the following. First, let us mention that 
the r = 1 and r = 2 parts in a free theory decouple and give, separately, identical results for observables due to the 
charge conjugation symmetry. They belong to two non communicating opposite charged mirror universes. Therefore 
the differences between the 8 component theory and the 4 component one can not show themselves if the r = 1 and 
r = 2 parts are not produced in coherent states of the form 

\<t>) =cos^|l) + e i<5 sin0|2) (15) 

with <fi ^ and S arbitrary, which does not separat apart during its evolution in time. The following is a scenario in 
which quantum interference effects of the 8 component theory that are not present in the 4 component one can be 
observed for interacting theories. 

The possibility of the existence of vacuum phases in a relativistic massless fermion system induced by a condensation 
of fermion pairs and antifermion pairs and their possible relevance to physical hadronic system are studied in Refs. 
T2| . Suppose that the /? and oj phases discussed in Rcfs. [| 12 that mix the r = 1 and r — 2 mirror universes are 



produced in, e.g. a heavy ion collision, or an astronomical object, then a quantum state of the form given by Eq. Ill 
with a fixed i5 and <j> ^ can be produced in the collision region Such a state contains r = 1 and r — 2 parts 
which have the same 3 momentum p. After this state leaves the collision region where the density is effectively zero 



to propagate to the detector, they remains in a mixed state of the form given by Eq. 15 with r — 1 and r = 2 parts 
having the same 3 momentum, which guarantees that these two parts can propagate together without been spatially 
separated on their way to the detector. To distinguish these two components of the particle in the detector, one can 
generate a non zero electric potential A by, say, putting the detector into a hollow charged conductor. Inside this 
A ^ region, the r = 1 and r = 2 pieces begin to have a different time dependence exp[— i(E p =p qA°)t] with q the 
charge of the particle, which can generates detectable beats in the counting rate. To show this more explicitly, let us 
consider the counting rate of finding a particle in a space volume Af2; Paq = Jaq d 3 x \ (x\p, </>) | 2 is 



d 3 " 



x 
An 



cos 2 0|(x|p,l)r + sin^|(x|p,2)r 

+ sin <j) cos (j) ((x | p,l)(x | p,2)*e~ i& +c.c)] , (16) 

where p is the 3 momentum of the excitation, c.c. denotes complex conjugation and the single particle physical state 
at x satisfies (a; \ P p hys — ( x \- It is ready to see that the third term depends on time with an universal beat frequency 
2qA° for all directly pruduced particles with the same q when <j> ^ 0. This effect is absent if the conventional 4 
component theory is used. 

The other excitation mode with r = 2 in Eq. ^ that corresponds to particle has a different spin structure and 
momentum range from the r = 1, e_ excitation mode. This excitation mode is not expected to has any interference 
effects with the r = 1, e_ excitation mode. Their effects on physical observables are addative on the probability level 
instead of the amplitude level. So they simply reproduce the effects of those in the 4 component theory due to charge 
conjugation symmetry of the system. 

In summary, we developed an 8 component theory for the relativistic finite density systems that has better short 
distance behavior than the corresponding 4 component one. The introduction of a mirror universe is found to be 
necessary. The link between the possible mirror universe and ours can be established and remains to be searched 
for in the real world phenomena of heavy ion collisions, mechanism for the CP violation, baryogenesis in the early 
Universe, the origin of dark matter, etc. 
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FIGURES 



FIG. 1. The set of contours belonging to the same FMS class. Contour Cfms is the original FMS contour in the Minkowski 
space. Contour Ce is the Euclidean contour. Contour C+ is the quasiparticle contour. Azin denote the imaginary part of the 
integrand along the edges of its cuts (thick lines) on the physical po plane. 



FIG. 2. The set of contours that give the thermodynamics. Contour Co encloses the Matsubara poles of the integrand. 
Contours C+ and C_ are the ones needed in a real time formulation of the thermal field theory in the 4 component representation 
for fermion fields. 
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